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Abstract 



We present a detailed quantitative discussion of the measurement of the leptonic mixing 
angle sin^ 2^^i3 with a future reactor neutrino oscillation experiment consisting of a near 
and far detector. We perform a thorough analysis of the impact of various systematical 
errors and compare the resulting physics potential to the one of planned first-generation 
superbeam experiments. Furthermore, we investigate the complementarity of both types of 
experiments. We find that, under realistic assumptions, a determination of sin^ 2^13 down 
to 10^^ is possible with reactor experiments. They are thus highly competitive to first- 
generation superbeams and may be able to test sin^ 26*13 on shorter timescales. In addition, 
we find that the combination of a KamLAND-size reactor experiment with one or two 
superbeams could substantially improve the ability to access the neutrino mass hierarchy 
or the leptonic CP phase. 
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1 Introduction 



There exists now strong evidence for atmospheric [1,2] and solar neutrino oscillations [3]. 
The recent KamLAND reactor experiment [4] together with all the solar neutrino data 
clearly confirms the LMA solution for the mass splittings and mixings and excludes other 
leading flavor transition mechanisms (for a recent review see, for example, Ref. [5]). Fur- 
thermore, the CHOOZ reactor experiment [6,7] currently provides the most stringent upper 
bound for the small sub-leading parameter sin^ 2^^i3. Since a non- vanishing value of sin^ 2^13 
is a prerequisite for genuine three-flavor effects, such as leptonic CP violation, it is now one 
of the most important challenges to establish sin^ 26'i3 > or at least to improve the sen- 
sitivity limit. Some improvement on sin^ 2^13 may be obtained from conventional beam 
experiments, such as the ongoing K2K experiment [8], or the MINOS [9] and CNGS [10] 
experiments under construction [11]. Even better limits can be obtained with planned su- 
perbeams [12-18] or future neutrino factory experiments (for a summary, see Ref. [19]). 
Especially superbeam experiments, such as the JHF to Super-Kamiokande [12] and the 
NuMI [13] experiments, are being planned with the main purpose to measure sin^ 2^13. 
Both superbeam and neutrino factory experiments suffer, however, from the presence of 
parameter correlations and degeneracies [17,20-22], and the sin^2^i3 measurements would 
not be as good as expected from statistics and systematics only [23]. Therefore, various 
suggestions have been made to resolve the correlations and degeneracies by the combination 
of at least two long-baseline experiments [14-16,24,25]. 

Oscillation experiments at nuclear power plants have recently been suggested as an inter- 
esting alternative to superbeams for the measurement of sin^ 2^13 [26-29] . Reactor experi- 
ments have a long history in neutrino physics. Starting from the legendary Cowan-Reines 
experiment [30], many measurements at nuclear power plants have provided valuable infor- 
mation about neutrinos. Very important are the results of the Gosgen [31], Bugey [32], Palo 
Verde [33], and CHOOZ [6] experiments, which have lead to stringent hmits on electron 
antineutrino disappearance. Reactor neutrino experiments have become very prominent 
again due to the outstanding results of the KamLAND experiment [4] . For a recent review 
on reactor neutrino experiments, see Ref. [34]. 

In order to improve the information on sin^ 2^^i3, it has been proposed to use a reactor experi- 
ment with a near detector very close to the reactor complex and a far detector at a distance 
of ~ 1 km — 2 km. Systematical errors can be reduced in this way, and a sensitivity to 
sin^ 2^13 down to ~ 0.02 might be reachable [27-29]. Moreover, such a measurement would 
not be spoilt by correlations and degeneracies. In this study, we analyze such reactor exper- 
iments and compare their physics potential to the one of first-generation superbeams. For 
that purpose, we thoroughly investigate the impact of various systematical errors. In addi- 
tion, performing a separate and combined analysis of reactor and superbeam experiments in 
the general three flavor framework, we study the competitiveness and the complementarity 
of the physics potentials of these types of experiments. 

The outline of the paper is as follows. In Section 2, we present the framework of neutrino 
oscillations and give a qualitative discussion of the sin^ 2^13 measurement at reactor and 
superbeam experiments using analytical formulas for the oscillation probabilities. In Sec- 
tion 3, we describe in detail how we simulate the experiments. In Section 4, we present 
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our results on the sensitivity limit on sin^ '2.9i^ obtainable at a reactor including a thorough 
discussion of various systematical uncertainties. In Section 5, we compare the limit as well 
as the accuracy for sin^ 26*13 obtainable at reactors with the one of superbeam experiments, 
whereas in Section 6 we show how a combined analysis of reactor and superbeam experi- 
ments can improve significantly the possibilities to identify the neutrino mass hierarchy and 
to discover leptonic CP violation. Our conclusions are presented in Section 7. Details of the 
statistical analysis are given in Appendix A, and in Appendix B we discuss experimental 
details of reactor experiments, including a summary of the key assumptions adopted in this 
work. In Appendix C we investigate the impact of the position of the near detector on the 
obtainable sensitivity. For readers mainly interested in our physics results, we recommend 
to proceed directly to Sections 4.1, 5, 6, and 7. 

2 The framework of neutrino oscillations 

Our results are based on a complete numerical three-flavor analysis. However, it is use- 
ful to have a qualitative analytical understanding of most effects. Therefore, we expand 
the relevant oscillation probabilities in terms of the small mass hierarchy parameter a = 
Am|]^/Am3i and the small mixing angle sin 26*13 using the standard parameterization of the 
leptonic mixing matrix U [35]. For the superbeams considered in this paper, we use, for 
the sake of simplicity, the expansion from Refs. [36-38] in vacuum, which is at least a good 
approximation for the JHF to Super-Kamiokande experiment.^ For the appearance signal 
with terms up to the second order, i.e., proportional to sin^ 26*13, sin 26*13 " ol-i and a^, one 
has 

P^e - sin^ 26*13 sin^ 6*23 sin^ A31 

^ a sin 2^13 sin 5cy> cos ^13 sin 2^12 sin 2^23 sin^ A31 

— a sin 2^13 cos t^cp cos ^13 sin 2^12 sin 2^23 cos A31 sin^ A31 

+ cos2^23sin2 2^i2sin2A3i, (1) 

where Ajj = AmfjL/{4E) = (mf — m'j) L / {4E) , and the sign of the second term refers to 
neutrinos (minus) or antineutrinos (plus). Depending on the actual values of a and sin 2^13, 
each of the individual terms in Equation (1) gets a relative weight. Thus, since within the 
LMA-allowed allowed region this relative weight can favor different terms in Equation (1), it 
turns out that the sin^ 2^i3-Am|^-plane is appropriate to illustrate the experimental poten- 
tial to measure sin^2^i3, the mass hierarchy, and CP effects. For example, for large a and 
large sin^ 2^^i3 the second and third terms are favored, which means that CP measurements 
become possible especially for large Am^i. On the other hand, for small a the first term 
in Equation (1) can be measured in a clean way without being affected by the other terms. 
Therefore, sin^ 2^13 or the sign of Awg^, entering the first term through matter effects, can 
be accessed there. In fact, it can be shown that a simultaneous measurement of Sqp and 
the sign of Amg^^ independent of the true value of Am^i is hardly possible for the first- 
generation superbeams or their combination [15]. The main reason for that problem is that 

^For similar expansions in matter, see, for example, Refs. [37-39]. 
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supcrbcams suffer from parameter correlations and degeneracies coming from tlie different 
combinations of parameters in Equation (1). Degeneracies are defined as solutions in pa- 
rameter space disconnected from the best-fit region at tlie chosen confidence level. Many 
of the degeneracy problems originate in the summation of the four terms in Equation (1) 
especially for large a, since changes of one parameter value can be often compensated by 
adjusting another one in a different term. This leads to the {6,613) [21], sgn{Amli) [17], 
and {023,7c/2 — 623) [20] degeneracies, i.e., and overall "eight-fold" degeneracy [22]. For 
superbeams, the {5, 6'i3)-degeneracy does no appear as a disconnected solution. In addi- 
tion, we choose the atmospheric best-fit value 623 — 7r/4, which means that especially the 
sgn(Am3i)-degeneracy will affect our results. We include the correlations and degeneracies, 
unless otherwise stated, in our results, and discuss their influence where appropriate. A 
detailed illustration of different sources of measurement errors and their impact can, for 
example, be found in Ref. [23]. In addition, the role of the degeneracies and the potential 
to resolve them has, for example, been studied in Refs. [14-16,24,25]. 

For the reactor experiments, we can, for short baselines, safely neglect matter effects. We 
find, up to second order in sin 2^13 and a. 



At the first atmospheric oscillation maximum, A31 is approximately 7r/2 and sin^Asi is 
close to one, which means that the second term on the right-hand side of this equation is 
also very small for sin^ 2^13 > 10^'^ and can for many purposes be neglected. The reactor 
measurement is at short baselines for large enough sin^ 2^^i3 therefore dominated by the 
product of sin^ 2^13 and sin^ A31, which must be measured as deviation from one. The 
simple structure of Equation (2) implies that, in comparison to the superbeams, correlations 
and degeneracies only play a minor role in reactor experiments - they are almost exclusively 
dominated by statistical and systematical errors. In other words, they can be used as "clean 
laboratories for ^13 measurements" [28]. Especially, the behavior in the sin^ 2^i3-Am2i-plane 
will be different to the one for superbeams, since Equation (2) is almost independent of 
Aml^. The direct comparison of reactor experiments and superbeams with respect to the 
most important dependencies will be one of the main aspects of this study. Equation (2) 
also demonstrates the limitations of reactor experiments, since there is no dependence on 
^23 1 ^CPi si'iid the sign of Amg-,^ in this formula. However, as we will show, the sensitivity 
to the sign of Am|j^ and to CP violation of superbeam experiments can be significantly 
improved by the combined analysis of reactor and superbeam experiments, through the 
precise determination of sin^ 2^13 at the reactor. 

All results within this study are, unless otherwise stated, calculated for the following values 
for the oscillation parameters, given with the currently allowed 3(T-ranges: 



1-P. 



ee 



~ sin^ 26*13 sin^ A31 + o? A31 cos^ 6*13 sin^ 26*12 ■ 



(2) 



Am\^ = 3^^10-^eV^ 

Sin2 2^23 = 1^0.2, 

Am^i = 7^f•10-^eV^ 

sin2 2^i2 = 0.81°:^ 



(3) 
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These numbers are motivated by recent global fits to atmospheric plus K2K data, such 
as in Refs. [2,8,40], and KamLAND plus solar neutrino experiments, such as in Ref. [41]. 
Throughout this work, we refer to the solar parameters given in Equation (3) as the LMA 
solution. In some cases, we will differentiate between the so-called LMA-I and LMA-II 
solutions with the best-fit values Am|i = 7-10'^ eV^ and Am|i = 1.4-10~^eV^, respectively, 
which both are covered by the 3cr-allowed range in Equation (3). Unless otherwise stated, 
we assume a normal mass hierarchy, i.e., Amg^^ = +3.0 ■ 10^^ cV^. For sin^2^^i3, we only 
allow values below the CHOOZ bound [6], i.e., sin^ 26*13 ^ 0.1. For the CP phase, we do 
not make specific assumptions, i.e., we allow any value between and 27r. 

3 The experiments and their analysis 

In order to reliably assess the physics potential of a given experiment, a realistic event 
rate calculation in connection with a proper statistical treatment of the simulated data 
is needed.^ The calculation of event rates is basically a convolution of the flux spectrum 
with the cross sections, the oscillation probability, the detector efficiency, and the detector 
energy response function. The resulting events rates are then the basis for a x^-analysis, 
where systematical uncertainties, correlations and degeneracies are properly included. The 
analysis for the superbeam experiments is done as in Refs. [15,23] . It is based on a Poissonian 
X^, since the appearance channel may have very low event numbers. On the other hand, 
the event rates in the disappearance channel of reactor experiments are quite large and we 
can use a GauBian x^, which has the advantage of allowing a more transparent inchision of 
systematical errors. For the detailed definition of the functions, we refer to Appendix A. 
Eventually, the combined analysis of reactor and superbeam experiments is done in a similar 
way as in Ref. [15], where we pay special attention to the correct treatment of degeneracies. 
All of the figures in this work are shown at the 90% confidence level. 

Once we have obtained a x^-valuc including the effects of systematical uncertainties, we 
project onto the parameter of interest in the six-dimensional space of the oscillation pa- 
rameters Amli, Am2i, 6*12, 6*13, 6*23, and 6cp- In addition, we take into account external 
information coming from other experiments, i.e., we assume to know the parameters Am|^, 
sin2 2^i2, and Ami^ with a 10% la-error from K2K [8], MINOS [9], CNGS [10], and Kam- 
LAND [4,42] results. This should be realistic at the time when the reactor or superbeam 
analysis will be performed. It turns out that our results are however rather insensitive 
towards the precise values of those errors as long as they are reasonably small, i.e., below 
about 30% to 50%. The matter density for the superbeam experiments is assumed to be 
known to 5% [43], which is essentially only relevant for the NuMI experiment. In the follow- 
ing two subsections, we provide further details on the simulations of the reactor experiments 
with near and far detectors and superbeams. 

^Further details about the general calculation and analysis techniques used in this work can be found in 
Ref. [23]. 
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3.1 Future reactor experiments with near and far detectors 

Nuclear fission reactors are a strong and pure source of low energy u^- Appearance exper- 
iments are not possible due to the low energies, and the inverse /3-decay with an energy 
threshold of 1.804 MeV is by far the dominant detection reaction: 

+ p —>■ e'^ + n . (4) 

This reaction has a very distinctive experimental signature which consists of the 7-rays 
from the annihilation of the positron and of the delayed signal of the neutron capture. This 
delayed coincidence allows to reject most of the background events. The relation between 
the positron energy and the neutrino energy is given by 

Ei^^ = Ee+ + {rrin - nip) + 0{Ej^Jmn) , (5) 

where Eg+ is the sum of the kinetic energy of the positron and the rest mass of the positron. 
The energy visible in the detector is determined by Ey-^g — Eg+ + 511 keV, where the ad- 
ditional 511 keV come from the mass of the electron which with the positron annihilates. 
Thus, a neutrino with the threshold energy produces already 2 ■ 511keV of visible energy. 
A precise measurement of the visible energy E^i^ yields according to Equation 5 a unique 
determination of the neutrino energy E^^ (for details of the energy reconstruction, see Ap- 
pendix B). The cross section for inverse /3-decay has approximately the form 

(T{Eg+) ~ Pe+Ee+ , (6) 

where r„ is the lifetime of a free neutron and / is the free neutron decay phase space 
factor. In our numerical calculations we use the cross sections from Ref. [44] including 
higher order corrections. For the neutrino flux spectrum from a nuclear power reactor we 
use the parameterization of Ref. [45] , and we adopt a fuel composition as in Ref. [4] . 

We assume a detector technology very similar to the CHOOZ [7] and KamLAND [4] de- 
tectors, as well as the Borexino counting test facility [46] (c/.. Appendix B). In this work, 

we suggest a reactor experiment designed to measure ^13. This measurement is based on 
the ability to detect small spectral distortions in the positron event rates due to neutrino 
oscillations. For a sensitivity to sin^ 26*13 of the order of 10~^ this requires a measurement ac- 
curacy of about 1%. The total experimental uncertainties are, however, considerably larger 
and come from various sources, such as reactor burn up effects, average neutrino yields, and 
cross section. In order to control these uncertainties, it is necessary to use a near detector, 
which accurately determines the energy dependence and total normalization of the neutrino 
flux [27,28]. The details of the near detector simulation are described in Appendices A and 
C. We furthermore do not include backgrounds since, as we discuss in detail in Appendix B, 
they can be suppressed to a negligible level [47] . 

We do not propose a specific reactor site in this study and for simplicity we assume a 
single reactor block. For a given site with several reactors effects like different distances 
of the cores to near and far detectors or on-off times of the individual reactors have to be 
included. However, one expects that the main results of this work will also apply in such 
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Figure 1: The sin^ 26*13 sensitivity limit as function of the baseline L for the Reactor-I (left 
plot) and Reactor-I I (right plot) setups. It is shown for several values for Amli at the 90% 
confidence level. The vertical lines mark the baseline we are using for our standard setups. 



more complicated situations. We define the integrated luminosity £ of a reactor experiment 
in units of fiducial detector mass [tons] x thermal reactor power [GW] x running time 
[years]. ^ In Figure 1, the dependence of the sin^ 26*13 sensitivity on the baseline of the far 
detector is shown for two values of the integrated luminosity 400 1 GW y (left-hand panel) 
and 8000 tGWy (right-hand panel), as well as several values of Amg^. It turns out that 
a baseline of 1.7 km performs reasonably well for all of the considered cases, which is in 
agreement with Ref. [28]. This choice is also best for the given uncertainty of Am|^ within 
the atmospheric allowed region, as it lies very close to the optimum for all values of Am|j^. 
We assume that the near detector is located at around 0.17 km in order to avoid oscillation 
effects. At those very short distances, the contribution of other power stations is far below 
1% of the total rate. As long their contribution is known better than to about 10%, they 
do therefore not contribute to the total error to more than 0.1% and we can safely ignore 
them. 

Throughout this work, we use two experimental benchmark setups labeled as Reactor-I and 
Reactor-I I, which are defined in Table 1. The two values of the integrated luminosity in 
the table reflect the range of possible experiments with detector masses of the order of 
lot [27,28] up to the order of > 100 1 [4,47]^. The near and far detectors are assumed to 
be identical (maybe apart from their size) in order to minimize the impact of systematical 
errors. We assume as our standard value for the uncertainty on the event normalization 

"^Notc that our definition of the integrated luminosity assumes a 100% detection efficiency and no dead- 
times. In addition, it includes the fiducial detector mass, not the total detector mass. Thus, for a specific 
setup, the losses due to these factors imply a rescaling of the luminosity. 

^Typical reactors have a thermal power ~ 2GW. Note, however, that reactor stations with a total 
thermal power of up to ~ 24 GW exist [28] . 
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Reactor- 1 


Reactor- II 


Integrated luminosity 
Unoscillated events 

Ccal 

Baseline 


400tGWy 
31493 

0.8% 

0.5% 
1.7 km 


8000 tGWy 
629 867 
0.8% 
0.5% 
1.7 km 



Table 1: The most important parameter values of the two reactor benchmark setups 
Reactor-I and Reactor-ll used in this work. For more information see Appendix B. 

Cnorm = 0.8% [28]. Tliis lias to be considered as an effective error, receiving contributions 
from individual uncertainties of the near and far detectors, as well as from over-all flux 
uncertainties (for details, see Section 4.2 and Appendix A). In addition, we include the 
effect of an energy cahbration error (Jcai = 0.5% [7,48] 

3.2 First-generation Superbeams 

The two superbeam experiments considered in this study are the JHF to Super-Kamiokande 
[12] and the NuMI [13] proposals, referred to as "JHF-SK" and "NuMI". Both projects will 
use the i/^ appearance channel. Neutrino beams which are produced by meson decays 
always contain irreducible fractions of z/g, z/g, and z/^j contaminations, as well as they have 
a large high energy tail. Both experiments use therefore an off-axis setup to make the 
spectrum much narrower in energy and to suppress the beam contaminations [49]. An off- 
axis beam reaches in this way the low background levels necessary for a good sensitivity 
to the Ve appearance signal. Both experiments are planned to be operated at nearly 

the same L/ E. which is optimized for the first maximum of the atmospheric oscillation for 
Am?i = 3.0- 10-3 eV^ 





JHF-SK 


NuMI 


Beam 


Baseline 


295 km 


712 km 


Target Power 


0.77 MW 


0.4 MW 


Off-axis angle 


2° 


0.72° 


Mean energy 


0.76 GeV 


2.22 GeV 


Mean L/E 


385kmGeV-^ 


320kmGeV-^ 


Detector 


Technology 


Water Chcrenkov 


Low-Z calorimeter 


Fiducial mass 


22.5 kt 


17 kt 


Running period 


5 years 


5 years 



Table 2: The two superbeams and their detectors as given in Refs. [12, 13] 

Due to the different energies of the two beams, different detector technologies are used. 
For the JHF beam, Super-Kamiokande, a water Cherenkov detector with a fiducial mass 
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of 22.5 kt, is assumed. The Super-Kamiokande detector has an excellent electron muon 
separation and NC (neutral current) rejection. For the NuMI beam, a low-Z calorimeter with 
a fiducial mass of 17kt is planned, because the hadronic fraction of the energy deposition 
is much larger at those energies. In spite of the very different detector technologies, their 
performances in terms of background levels and efficiencies are rather similar. The actual 
numbers for these quantities and the respective energy resolution of the detectors can be 
found in [15]. 

In Ref. [15], several modifications of those setups have been considered in order to improve 
either the sensitivity to CP violation or to the mass hierarchy. There, it is demonstrated 
that JHF-SK with about 1/4 of the total running time with neutrinos and about 3/4 with 
antineutrinos performs best for the CP measurement^, a setup which is labeled as JHF-SKcc 
in this work. However, for the determination of the mass hierarchy the NuMI setup with an 
increased baseline of 890 km at an off-axis angle of 0.72° performs better because of larger 
matter effects, which will be labeled as NuMI@890km in this work. 

4 Measuring sin^ 2^13 at a reactor, and the impact of systematical 
errors 

In this section, we investigate in detail the potential of a reactor experiment to measure 
sin^26'i3. In Section 4.1, we discuss obtainable limits on sin^26'i3, and we present in Sec- 
tions 4.2 and 4.3 a detailed discussion of systematical errors. The results of these sections are 
obtained by neglecting correlations of the oscillation parameters, i.e., all parameters except 
from sin^ 2^13 are fixed to the values given in Equation (3). However, as it can be inferred 
from the discussion related to Equation (2) and as we will see explicitly in the numerical 
calculations presented in Sections 5 and 6, the infiuence of other oscillation parameters on 
the sin^ 2^13 sensitivity limit from reactors is very small. 

4.1 The sensitivity limit to sin^2^i3 

We define the sensitivity or sensitivity limit to sin^ 26'i3 as the largest value of sin^ 26*13, 
which is allowed at a given confidence level if the true value is sin^ 2^13 = 0. In Figure 2, 
we show the sensitivity to sin^ 2^13 as a function of the integrated luminosity C in units 
of detector mass [tons] x thermal reactor power [CW] x running time [years]. In this 
figure, the lower diagonal curve corresponds to the idealized case of statistical errors only, 
and shows just the expected l/\/~C scaling. Our standard values for the errors on the 
normalization anorm = 0.8% and energy calibration Ucai = 0.5% lead to the thick curve, which 
represents one of the main results of this paper. At a luminosity around lOOtGWy, we 
detect a departure from the statistics dominated regime into a fiatter systematics dominated 
region. As discussed in Sections 4.2 and 4.3, this effect is dominated by the error on 
the normalization Unorm, whereas the energy calibration error (Jcai only plays a minor role. 
However, at large luminosities > lO^tGWy, the slope of the curve changes, and we are 
entering again a statistics dominated region with a 1/ y/C scaling. This interesting behavior 

^This setup has about equal total numbers of neutrino and antineutrino events. 
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Figure 2: The sensitivity to sin^ 2^^i3 as a function of the integrated luminosity for different 
values of the normalization error cXnorm o^nd the energy calibration error cXcai at the 90% 
confidence level, where the thick curve corresponds to our standard values. The vertical 
lines mark the luminosities of the CHOOZ experiment and our standard setups Reactor- 1 
and Reactor-ll as defined in Table 1. The horizontal line shows a typical sin^ 2^13 sensitivity 
limit obtainable by the JHF-SK superbeam experiment for the same parameter values (taken 
from Section 5). 

can be understood as follows: In principle, the far detector measures some combination of the 
flux normalization and sin^ 2^13. The turnover of the sensitivity line into the second statistics 
dominated region occurs at the point, where the determination of the normalization by the 
far detector itself becomes more accurate than the external input o"norm- This means that 
the limit becomes insensitive to the actual value of o"norm- We illustrate this by the upper 
thin black line, which shows the luminosity scaling for the case of larger systematical errors. 
As an example, we choose values of iTnorm = 1.7% for the normalization and (Jcai = 1.1% 
for the energy calibration. We find that, in this case, the transition to the systematics 
dominated regime occurs at much smaller luminosities. However, for large luminosities, the 
same limit is approached as for the more optimistic case. The diagonal gray curve shows 
the limit for no constraint at all on the normalization and energy calibration.® Even in this 
extreme case, we obtain the same limit for high luminosities. 

This discussion demonstrates that the magnitude of the systematical error determines the 
position of the sensitivity plateau, but essentially does not affect the sensitivity at large 
luminosities. We hence conclude that, for the case of the Reactor-I setup, the systematical 

^Although we leave the normalization free in the fit, we assume that the shape is known. 
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normalization error dominates. In order to obtain a reliable limit, it should therefore be well 
under control. For large luminosities, such as for the Reactor-ll setup, the sensitivity limit 
is independent of normalization errors because of the second statistics dominated regime. 
The robustness of this result will be further discussed in Section 4.2. 

The horizontal line in Figure 2 shows the typical sin^26'i3 sensitivity limit, which can be 
obtained from the first-generation JHF-SK superbeam experiment including correlations and 
degeneracies. It can be inferred that even the Reactor- 1 experiment gives comparable hmits. 
The competitiveness and complementarity of the information from reactor and superbeam 
experiments will be discussed in greater detail in Sections 5 and 6. Finally, we note that 
with a KamLAND-like detector of 1 kt at a 10 GW nuclear reactor with a running time 
of ~ 5 years, the high statistic region of a luminosity between lO^tGWy and lO^tGWy 
seems to be reachable. 

4.2 The impact of systematical errors 

We discuss now the effects of various systematical errors on the obtainable sensitivity to 
sin^ 2^13, where technical details are given in Appendix A. Since it is a difficult task to 
estimate systematical errors of a future experiment, our strategy has been to find realistic 
values for various normalization errors {e.g., total neutrino fiux uncertainty, fiducial mass 
uncertainty) and energy calibration errors. The chosen numbers are guided by the values 
obtained in existing reactor experiments, such as the CHOOZ [6, 7] or KamLAND [4] ex- 
periments. For other types of errors, such as the shape uncertainty or the experimental 
systematical error, we have adopted a conservative approach by choosing the worst case 
situation of completely uncorrelated errors. In detail, we consider the following effects (for 
the exact definitions, see Appendix A): 

1. We take into account a common overall normalization error for the event rates of the 
near and far detectors. Such an error could, for example, come from the uncertainty of 
the neutrino fiux normalization or the error on the detection cross section. Typically, 
it is of the order of a few percent. 

2. We include uncorrelated normalization uncertainties of the near and far detectors. 
Here contributes, for instance, the error on the fiducial mass of each detector. We 
assume that in this case an error below 1% can be reached. 

3. We take into account the energy calibration uncertainty by introducing a parameter 

for each detector {A = N, F), and replace the observed energy i?obs by {l + g'^)Eobs- 
We assume that the energy calibration is known within an error of cTcai ~ 0.5%. 

4. In order to take into account an uncertainty of the shape of the expected energy 
spectrum, we introduce an error o"shapc on the theoretical prediction for each energy 
bin which is completely uncorrelated between different energy bins. This corresponds 
to the most pessimistic assumption of no knowledge of possible shape distortions. 
However, we choose this error fully correlated between the corresponding bins in the 
near and far detector, since shape distortions should affect the signals in both detectors 
of equal technology in the same way. 
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5. We include the possibility of an uncorrelated experimental systematical error a^xp- 
Such an error could, for example, result from insufficient knowledge of some source 
of background. We call this uncertainty "bin-to-bin error" and take it completely 
uncorrelated between energy bins, as well as between the near and far detectors. Note 
that this corresponds again to the worst case scenario, and values of (Jexp at the per 
mill level should be realistic. 

As it is explicitely demonstrated in Appendix A, the overall normalization error, the in- 
dividual normalization errors of the two detectors, and the energy calibration error of the 
near detector can be merged into an effective normalization error Unorm for the far detector. 
Assuming realistic vahics of 2% [7] for the total normalization uncertainty and 0.6% [48] for 
the detector-specific uncertainty, we obtain from Equations (18) and (19) in Appendix A an 
effective normalization error of (Tnorm = 0.8%, which is the standard value for our numerical 
calculations. 

Luminosity Scaling 

0.03 II I I I I I I I I I I I I I I I I I I I III 




10^ 10^ lO'* 10^ 

Integrated luminosity £ [t-GW y] 



Figure 3: The sensitivity to sin^ 26'i3 o,s a function of the integrated luminosity for several 
selected error sources at the 90% confidence level. The thick curve corresponds to our stan- 
dard values for the errors on the flux normalization and the energy calibration. The gray 
curve has been calculated by adding an error of 2% for the shape of the predicted energy 
spectrum. For the curves labeled with "bin-to-bin" , we include an uncorrelated experimental 
systematical error as given by the corresponding label. The vertical lines mark the integrated 
luminosities of our standard setups Reactor-I and Reactor-ll as defined in Table 1. 

The impact of a spectral shape uncertainty and of experimental systematical errors is illus- 
trated in Figure 3, where we show the luminosity scaling of the sin^ 2^^i3-limit for some values 
of these errors. The curve labeled with "2% shape" corresponds to an uncorrelated shape 
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error on the predicted energy spectrum, such as it is described in item 4 above. This error 
covers a large class of systematical effects, such as the uncertainty on the /3-decay spectrum 
of various isotopes in the reactor, uncertainties of the fuel composition, or burn-up effects. 
As it can be inferred from the figure, such an error has very little impact on the sensitivity 
limit. The reason for this is that shape uncertainties can be reduced very efficiently by the 
near detector, assTiming that there are at least 10 times more events in the near detector 
than in the far detector. Our assumption of an uncorrelated error is the most pessimistic 
theoretical error one could imagine. In more realistic cases, one may expect some correlation 
of the shape uncertainty between the energy bins, which could be even better eliminated by 
the near detector. 

The presence of a systematical, experimental bin-to-bin uncorrelated error, as described in 
item 5, is potentially more problematic for the measurement than any other error, since it 
cannot be eliminated with the help of the near detector. As it is shown in Figure 3, such 
an error with a value below 0.5% has no effect for the Reactor-I setup, but it spoils the 
statistics dominated region at high luminosities and would be important for large detectors, 
such as Reactor-I I. However, we note that 0.5% for an uncorrelated error of this type is a 
relatively large number, which should be reducible down to the 0.1% level. In this region, 
the impact on the sin^ 2^13-limit becomes very small. Suppose, for example, that the number 
of background events in the detector is 1% of the reactor neutrino events. In this case, a 
10% knowledge of this background is sufficient to reach a 0.1% error on the experimental 
event rate. 

4.3 Normalization and energy calibration errors, and the impact of an energy 
cut 

Since we have already demonstrated that the impact of an error on the expected energy 
shape is very small, we will neglect it for the rest of this work. Moreover, we further on 
assume that the experimental bin-to-bin uncorrelated systematical error can be reduced 
to the 0.1% level and hence can also be neglected. Then, as shown in Appendix A, all 
systematical effects considered in this paper can be reduced to two systematical errors, 
which are the effective normalization uncertainty cXnorm and the energy calibration error Ucai 
in the far detector. We will discuss now in this subsection the impact of these two errors 
in greater detail, combined with a possible low energy cut which could be required in some 
cases to eliminate backgrounds. 

In Figure 4, we illustrate the sin^ 2^13 sensitivity limits for various values of o"norm and 
(Jcai- The vertical lines at the left edge of each bar correspond to the idealized case of 
statistical errors only. The bars indicate, how the limit deteriorates by including either the 
normalization error or the calibration error or both. Let us first focus on the cases labeled 
"no cut", where the full energy range of 1.8 MeV < E^, < 8MeV is considered. The final 
sensitivity limit is dominated by the normalization error, and the calibration error has a 
rather small impact, although it is not negligible. We also show in Figure 4 the extreme case, 
where (Xnorm and cXcai are set to infinity, i.e., the event normalization and the energy scale 
are treated as free parameters in the fit (c/.. Equation (17) in Appendix A). We learn from 
this analysis that for (Tcai > only, even if we include the energy scale as a free parameter 
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Figure 4: The sensitivity to sin^ 2^13 for the Reactor-I and Reactor-I I setups at the 90% 
confidence level. The left vertical line in each bar corresponds to statistical errors only, 
whereas the right edge shows the limit for different combinations of normalization and cal- 
ibration errors. The sensitivity limits are shown for o"norm = 0.8% and cxcai = 0.5%, or for 
Cnorm = o"cai oo whcrc indicated explicitly. Furthermore, the sensitivity limits are shown 
without cut on the energy, and for an energy cut at > 3.4 MeV. 

in tlie fit, tiie limit is essentially the same as the pure statistics limit. This demonstrates 
again that an energy scale uncertainty is not very important. In addition, we observe that 
for Reactor-I the final sensitivity limit is significantly worse for (Tnorm = CTcai ^ 00 than 
for cTnorm = 0.8%, (Teal = 0.5%, whcrcas the limits for Reactor-I I are essentially unchanged. 
This again illustrates the effect already discussed earlier in the context of Figure 2: At 
high luminosities the oscillation signal is dominated by the spectral information and the 
far detector itself provides an accurate measurement of the overall fiux, which means that 
the sin^ 2^13 sensitivity limit becomes insensitive to the value of Unorm- On the contrary, for 
lower luminosities, such as for the Reactor-I setup, a precise knowledge on (Tnorm is rather 
important, and the limit can be significantly improved for small normalization errors. 

The cases in Figure 4 labeled "with cut" correspond to the situation when a low energy 
cut is imposed. Such a cut could be motivated by some low energy background, such 
as from radioactivity of the detector material [7]. For illustration, we have chosen the 
cut > 3.4 MeV corresponding to i^vis > 2.6 MeV, as it is applied in the KamLAND 
experiment to avoid the contribution of geo-neutrinos. Apart from the slight worsening of 
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Figure 5: The ratio of the event rates at the 3(T limit of sin^ 2^^i3 to the event rates for 
sin2 2^i3 = for Reactor-ll (sin^ 2^13 = 1.5 ■ 10"^ in the left plot and s\-n^2ei^ = 2.5 • 10"^ 
in the right plot). In the left-hand panel, the full energy range is included, whereas in 
the right-hand panel, an energy cut at E^j > 3.4 MeV is imposed. The thin solid curves 
correspond to the case of no systematical errors, whereas the long-dashed curves correspond 
to a calibration error of cXcai = 0.5%, the short-dashed curves to a normalization error of 
Cnorm = 0.8%, and the thick solid curves to both errors switched on simultaneously. The 
shaded regions show the statistical errors, and the numbers given for each curve correspond 
to the separate contribution to the total value of Ax^ of the data and the systematical errors, 
see Equation (17). 

the statistics-only limit due to smaller event numbers, we find that the cut has a rather 
drastic impact for the large detector Reactor-ll. The reason for this effect is illustrated in 
Figure 5, where we show the normalized energy spectrum for the following four cases: no 
systematical errors, either (Jnorm or Ucai included, and both iTnorm and Ucai simultaneously 
included. One can see that the cut at 3.4 MeV is close to the oscillation minimum at this 
baseline of 1.7 km. For such an unfortunate choice, the interplay of normalization and 
calibration can very efficiently reduce the oscillation pattern in the spectrum by making 
the signal essentially flat. In fact, we have found that, in the presence of an energy cut, 
our standard values for (Jnorm and (Tcai are sufficient to completely destroy the statistics 
dominated regime at high luminosities, similar to the "bin-to-bin" error shown in Figure 3. 
This problem can be avoided by making sure that the oscillation minimum is safely included 
inside the accessible energy range, such as for the case when no cut is applied. With this 
choice, it is impossible to destroy the oscillation signature by shifting the normalization or 
stretching the energy scale, as it is clear from the left plot of Figure 5. If a low energy cut 
cannot be avoided because of backgrounds, one should eventually consider to change the 
baseline in order to make sure that the oscillation minimum is covered by the remaining 
energy window. 



14 



5 The measurement of sin^ 2^13: reactors versus superbeams 

In this section, we compare the potential to mcasiu^c sin^ 2^13 in reactor experiments with 
the one of superbeams. As far as sin^ 29i^ is concerned, we will demonstrate that even the 
small reactor setup Reactor-I can compete with a first-generation superbeam, such as JHF-SK 
or NuMI. There are essentially two interesting measurements for sin^2^i3: the sensitivity 
to sin^20i3 and the precision of the measurement of sin^2^i3. The discussion is therefore 
divided into two parts, where we also define the respective quantities. 

5.1 The sensitivity to sin^26'i3 

We define the sensitivity or sensitivity limit to sin^2^i3 as the largest value of sin^2^i3, 
which fits the true value sin^ 2^13 = at the chosen confidence level. For an experiment 
or a combination of experiments, it reflects the range between the sensitivity limit and 
the CHOOZ bound, where sin^ 26*13 > could be detected. With this definition, it should 
be clear that especially for superbeams the sgn(Am|i)-degeneracy has to be taken into 
account in the calculations, since the degenerate solution at — |Am|i| may allow a larger 
value of sin^ 2^^i3 fitting sin^ 2^13 = than the best-fit solution. In addition, it has been 
demonstrated in Ref. [15] that with this definition the sensitivity limits for the normal and 
inverted mass hierarchies are equal for the superbeams, since the zero rate vectors for the 
appearance channel are equal for sin^ 2^13 = 0. Similarly, Equation (2) does not depend on 
the sign of Am|]^, which means that the sensitivity limit does not depend on the type of the 
hierarchy in this case, either. 

Figure 6 shows the sensitivity limits for JHF-SK as example for the superbeams and Reactor-I, 
as well as their combination at the 90% confidence level. In this figure, the left plot repre- 
sents the Amg^-dependence and the right plot the Amgi-dependence. In order to compare 
experiments with similar capabilities, we have chosen the Reactor-I setups for these plots. 
The Reactor-I I setup would be much better and would dominate the result. 

For the best-fit parameters used in this work, we obtain a sensitivity hmit of sin^ 2^13 ~ 
1.8 • 10~^ for Reactor-I at the 90% confidence level, as it can be read off from the figure. 
This result is in rather good agreement with the setups in Refs. [28] and [27], where similar 
values for the systematics uncertainties were used. In Figure 6, the Amg^^-dependence clearly 
refiects the fact that both the JHF-SK and Reactor-I experiments are optimized for a value 
of Ato|i ~ 3.0 • 10-^ eVl However, as we have seen in Figure 1, the broad reactor spectrum 
does not make this optimization peak as sharp as in the superbeam case, where the off- 
axis technology is actually used to produce a narrow-band beam. The risk of choosing 
a non-optimal value for the Am3^ -optimization within the atmospheric allowed region is 
therefore much lower in the reactor case, though the superbeam is marginally better at the 
best-fit value. From the combination of the two experiments, we do not observe synergy 
effects which would improve the performance beyond a simple addition of statistics. As it 
is demonstrated in Figure 3 of Ref. [15], the NuMI setup would lead to very similar results. 

As far as the Am^i-dependence is concerned, superbeams strongly suffer from the correlation 
with the CP phase. This can already be seen in Equation (1), where especially for large 
values of the hierarchy parameter a = Am|i/Am|i the second and third terms become 
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Figure 6: The sensitivity limits to sin^ 26i^ as functions of the true values of Am|]^ (left 
plot) or Arri^i (right plot), respectively. They are shown for the JHF-SK and Reactor-I setups 
as well as their combination at the 90% confidence level. In the left figure, the atmospheric 
excluded region is shaded in gray, and in the right figure, only the KamLAND-allowed region 
is shown. In both figures, the best-fit values are marked by the gray lines. 

larger, and sin^ 2^^i3 becomes highly correlated with 5cp- Thus, the JHF-SK setup loses 
almost any sensitivity to sin^ 2^13 below the CHOOZ bound for large values of /S.m\^ within 
the LMA-allowed region, as it can be seen in the right plot of Figure 6. On the other side, 
the reactor setup is hardly influenced by a, since for the short baseline the effects of a 
within the LMA-allowed region are small. Thus, though a httle bit worse for small values 
of Am|]^, Reactor-I has a much better sensitivity than JHF-SK in most of the LMA-allowed 
region. The combination of the two experiments does, as in the case of Amg^, not show 
significant synergy effects beyond adding the statistics of the two experiments. 

Apart from sin^2^^i3, the reactor experiment is sensitive to some parameters which can 
be more precisely measured by earlier experiments. We have therefore imposed external 
information on the parameters Am^^, Am2i, and sin^2^i2 appearing in Equation (2), and 
have studied the dependencies on the precision of this information. We find that the reactor 
experiment itself can measure the relevant parameters with a sufficient precision, as long 
as Am|^ is known better than to 30% — 50% from external measurements. This required 
external precision of Am|^ is easily achievable by conventional beam experiments, such 
as CNGS, MINOS, or K2K. The Reactor-I setup is in summary very competitive to the 
first-generation superbeams, especially since the risk of the unknown parameter values of 
Amg]^ and Amgi is considerably lower. With a higher luminosity, the Reactor-I I setup could 
moreover do much better than the first-generation superbeams on comparable timescales. 
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Figure 7: The relative error on log(sin^ 2^13) defined in Equation (7) as function of the 
true value o/sin^2^i3 /or JHF-SK, Reactor-I, Reactor-ll, and the combination o/JHF-SK and 
Reactor-I at the 90% confidence level. For experiments or combinations involving JHF-SK, 
the (by then unknown) true value of Sep leads to bands representing all possible values of 
6cp- In addition, for each value of 5cp, the worst among the best-fit and degenerate solutions 
was chosen. 

5.2 The precision of sin^ 2^13 

As soon as sin^ 26*13 7^ is established by an experiment, the obtainable precision of sin^ 26*13 
becomes interesting. This precision of the measurement of sin^ 2^13 depends, of course, on 
the true value of sin^ 2^13 itself. A meaningful quantity to discuss this precision is the 
relative error on log(sin^ 2^13) as a function of sin^ 2^13 [23], i.e., 

, / 2 . ^ rr..! Ilogfsin^ 2^^.^) - logfsin^ 2^S) I 

Rel. error on log sin^ 2^13) % = ^—^ 9 \ ^ • 100. 7) 

log(sin 26*13) 

Here logfsin 2efl) and log(sin2 29^1^) refer to the upper and lower intersections, respectively, 
of the x^-function with the 90% confidence level. With this definition, it should be obvious 
that this relative error becomes larger than 100% close to the sensitivity limit. It is shown in 
Figure 7 for JHF-SK and Reactor-I, their combination, and Reactor-I I. For any combination 
with a superbeam, we included in this calculation the degenerate solutions by taking the 
largest relative error of all degenerate solutions. In addition, we made no special assumptions 
about the unknown CP phase, which leads to the bands for experiments involving the 
superbeam in Figure 7. The lower edge of these bands corresponds to the best case, the 
upper edge to the worst case. This form of visualization takes into account that we do not 
know by then, which true value of 6cp has been realized by nature. The different results 
for different values of 5cp originate in the different shapes of the 5cp-^i3-correlation, which 
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means that this dependence does not necessarily imply that the experiment can measure 
Sep- 

From Figure 7, we find that the precision of sin^ 2^13 is much better for the reactor experi- 
ments for large values of sin^2^^i3, but for small values of sin^ 2^^i3 the supcrbeams become 
better. This behavior stems from the different nature of the signals: the reactor experiments 
measure ^13 in the disappearance channel, whereas the superbeams measure the appearance 
of events. For reactor experiments, the statistical and systematical errors are therefore basi- 
cally independent of sin^ 2^^i3, which implies that the precision of sin^ 2^13 vanishes for small 
values of sin^2^^i3. On the other hand, superbeams have, even for large values of sin^2^i3, 
only a limited number of signal events [23]. Therefore, the statistical errors are relatively 
large compared to a reactor experiment. For smaller values of sin^2^i3, however, the rel- 
ative statistical error decreases as 1/y/N. This results in a much better accuracy at low 
values of sin^2^^i3. The same loss of precision for reactor experiments can also be observed 
in the sin^ 2^^i3-Am|^-plane, as shown in Figure 2 of Ref. [28]. This figure illustrates that 
under a certain threshold value for sin^ 26*13 the shape of the sin^ 26'i3-Am|i allowed region 
changes dramatically. Another consequence of the shape of this region is that an external 
measurement of Amg^ does not influence the precision of sin^ 2^13 in Figure 7 very much, 
because for large values of sin^ 2^13 the parameters Am^^ and sin^ 2^13 are uncorrelated, 
and for small values of sin^ 2^^i3 the relative error as defined in this section quickly becomes 
very large before Amli and sin^ 2^^i3 become highly correlated. 

As far as the combination of Reactor- 1 and JHF-SK is concerned, there is some synergy for 
small values of sin^ 2^13. Though the precision of the reactor experiment deteriorates, both 
experiments together are somewhat better than the superbeam alone. The Reactor- 1 1 setup 
is however much better, even better than the combination of Reactor-I with JHF-SK for large 
values of sin^ 2^^i3, which means that a larger detector can really help in this case. After all, 
the superbeam suffers from the same problem as for the sensitivity limit: the larger Arn^i 
is, the worse the precision of sin^ 2^^i3 becomes. An unfortunate true value of the CP phase 
would additional erode the performance. Contrary to that, the Reactor-I I setup does not 
have these problems, which means that a reactor experiment with a large detector would be 
the optimal choice to measure sin^ 2^13 with lower risks coming from oscillation parameter 
uncertainties. 

6 The complementarity of reactor experiments and superbeams 

We have already seen that the reactor experiments are very competitive in measuring 
sin^ 2^13. There are, however, other parameters, which the reactor experiments cannot 
access satisfactorily, such as the sign of Am|i and CP violation. We will therefore demon- 
strate now, how the reactor experiments would fit in the larger picture of the next-generation 
experiments.^ We will use the Reactor-I I setup to identify the optimal combinations with the 
first-generation superbeams. This will demonstrate what reactor experiments can contribute 
in the best case. 

qualitative discussion of many interesting issues for the complementarity of future reactor and long- 
baseline experiments can also be found in Ref. [28]. 
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Figure 8: The relative error on Am|i in percent as function of the true value of sin^ 2^13 
/or JHF-SK, NuMI, Reactor-I, and Reactor-II. For the calculation of the error, the full width 
of the -function in the best-fit region is taken at the 90% confidence level. 

6.1 Precision measurements of the leading atmospheric parameters 

Among the interesting parameters for superbeams are the leading atmospheric parameters 
Am|]^ and sin^2^23- Compared to conventional beams, such as K2K, CNGS, and MINOS, 
the superbeams could achieve very high precisions for the leading atmospheric parameters. 
Reactor experiments are, at the short baselines we are considering, besides sin^2^^i3, only 
sensitive to Amli- This can be understood in terms of Equation (2): the atmospheric 
mixing angle 623 does not even appear in this equation and the sensitivity to the lead- 
ing solar parameters can only be sensibly achieved at longer baselines, where A21 is large 
enough. Therefore, the first-generation superbeams would supply a measurement without 
serious competition at least for sin^2^^23- The atmospheric mass squared difference Aml^, 
however, could be measured by superbeams as well as reactor experiments, as it is illus- 
trated in Figure 8. There is nevertheless one important difference between those two types 
of experiments. The superbeams measure Am^^ with the disappearance channels, which 
are dominated by the leading atmospheric oscillation, and the measurement thus hardly 
depends on the true value of sin^2^^i3. For the reactor experiments, though, the parameter 
Amg^ is part of the signal proportional to sin^ 26*13. Therefore, they are strongly affected 
by the true value of sin^ 2^13, as it is demonstrated in Figure 8. We hence conclude that 
long-baseline experiments are very important for precision measurements of the leading 
atmospheric parameters independent of the true value of sin^ 2^13. 
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Figure 9: The sensitivity to a positive sign of Arn^i (normal mass hierarchy) as function 
of the true values o/sin^2^i3 and Am2i within the KamLAND-allowed region. Sensitivity at 
the 90% confidence level exists on the right-hand sides of the curves. It is shown for several 
combinations of experiments: /or JHF-SK combined with NuMI at a baseline o/ 890 km [15], 
for NuMI at a baseline o/ 890 km combined with Reactor-ll, and for JHF-SK combined with 
NuMI at a baseline o/ 890 km and Reactor-ll. The LMA best-fit values are marked by the 
horizontal gray lines. 



6.2 The sensitivity to the sign of Am| 



As it has been demonstrated in Refs. [15,23] and elsewhere, the mass hierarchy is, due to 
the sgn(Am|^)-degeneracy, one of the hardest parameters to access for future long-basehne 
experiments. As it can be seen in Equation (1), the opposite sign of Aml^ opens, especially 
for large Am^i, the possibility of a degenerate solution at a different value of ^cpj which 
makes it very hard to determine the mass hierarchy. We define that an experiment is 
sensitive to a certain sign of Am|^, if there is no possible solution with the opposite sign of 
Am|^ below the chosen confidence level. This sensitivity depends, similarly to the precision 
of sin^ 26*13, on the unknown true value of 6cp- We take therefore the most conservative value 
of 5cp in order to show where sensitivity to the tested hierarchy exists independent of Sep- 
With this definition, it can be shown that neither the first-generation superbeams, such as 
NuMI or JHF-SK, nor their combination, do have any sensitivity to the mass hierarchy [15]. 
However, because of matter effects, which increase with the baseline, and the ability to 
resolve degeneracies with the combination of two superbeams [15,25], the combination of 
NuMI at a longer basehne together with JHF-SK is sensitive to the sign of Aml^ in a 
large region of the sin^ 2^i3-Am|i-plane. The result for a positive sign of Amg^ and the 
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combination of JHF-SK with NuMI at a baseline of 890 km^ from Rcf. [15] is, for example, 
shown in Figure 9. As it can be seen in the figure, the sensitivity to the normal mass 
hierarchy vanishes for large values of Am^i, which mainly comes from the sgn(Am|]^)- 
degeneracy. 

We find however that the combination of the Reactor-ll setup at the short baseline of 1.7 km 
with two superbeam experiments significantly improves the sensitivity to the mass hierarchy. 
Figure 9 demonstrates that this combination (thick solid curve) has a very good sensitiv- 
ity to Amg^ > which does not depend on ^m\^. The reactor experiment helps here 
indirectly to resolve the sgn(Am|^)-degeneracy by the precision measurement of sin^2^^i3. 
The combination of the three experiments thus leads to a much better sensitivity to the 
normal mass hierarchy at large values of Am2i within the LMA-allowed region. In order to 
demonstrate that one superbeam plus the reactor experiment at the short baseline is not 
enough to access the mass hierarchy, we show, in addition, in Figure 9 the combination of 
the superbeam with the larger matter effects, z.e., NuMI at 890km, together with Reactor-ll. 
Figure 9 is for a normal mass hierarchy, but the inverted mass hierarchy produces rather 
similar results [15]. 

In Ref. [50], it has been pointed out that for the case of a relatively large Am^i, a reactor 
neutrino experiment at a baseline of the order of 30 km might have some sensitivity to the 
sign of ^m\^ because of an interference term between Am^i and Amg^ (see also Rcf. [47]). 
This interference term is of third order proportional to a-sin^ 2^^i3. which means that it does 
not appear in our Equation (2). According to our numerical results, a determination of the 
mass hierarchy at the 90% confidence level based on this effect is not possible, even for a 
combination of two big reactor experiments with luminosities of lO^tGWy at baselines of 
1.7 km and 30 km. In order to observe this tiny effect, the ratio ^m^i/ ^m\-y and sin^ 2^13 
should be as large as possible, and the solar mixing has to be far from maximal mixing. 
Moreover, all these parameters have to been know with (unrealistically) high precisions. 

6.3 The sensitivity to CP violation 

There are various approaches to evaluate the capability of an experiment to access 5cp- One 
of those, which is quite intuitive to understand and qualitatively representative, is the sensi- 
tivity to maximal CP violation 5cp = ±7r/2. In this spirit, we define that an experiment (or 
a combination) is sensitive to maximal CP violation, if the true value 5cp = ±7r/2 does not 
fit the CP conserving values and tt at the chosen confidence level. With this definition, it 
is obvious that also the degenerate solutions have to be tested, since any degenerate solution 
with (5cp = or TT fitting the maximal CP violation would destroy the sensitivity. It has 
been demonstrated in Ref. [15] that for the first-generation superbeams the results quali- 
tatively do not depend very much on the choice of the normal or inverted mass hierarchy 
and the true value 5cp — +7r/2 or 5cp = — 7r/2. We show therefore only the results for 
the normal hierarchy and 5cp = +7r/2 in this work. In addition, it has been discussed in 
Refs. [15,25] and elsewhere that for the detection of CP violation with superbeams a com- 
bined neutrino and antineutrino running is very important - cither within one superbeam 

^This NuMI baseline corresponds to the longest allowed baseline for the already fixed decay pipe at an 
off-axis angle of 0.72°. 
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Figure 10: The sensitivity to maximal CP violation Sqp — +7r/2 as function of the true 

values ci/sin^2^i3 and Aml;^ within the KamL AND- allowed region. The sensitivity is at the 
90% confidence level given on the upper sides of the curves. It is shown for SV\¥SV\ (neutrino 
running only), JHF-SKcc (same overall running time split into neutrino and antineutrino 
running with about equal numbers of events), and JHF-SK (neutrino running only) combined 
with Reactor-ll. The LMA-II best-fit value is marked by the horizontal gray line. 

experiment or in the combination of different superbeams with different polarities. This 
result is not very surprising, since those two channels have a complementary dependence 
on the CP phase (c/., Equation (1)) and their combination helps to resolve the correlation 
between sin^2^i3 and 5cp- The antineutrino running has, however, one disadvantage: in 
order to accumulate similar event rates for comparable statistical weights, the antineutrino 
mode has to be operated much longer than the neutrino mode to compensate the lower cross 
section of antineutrinos. As illustrated for different setups in Figure 10, there is a possibil- 
ity to circumvent this problem with JHF-SK. Instead of resolving the correlation between 
sin^ 26*13 and (5cp with the antineutrino channel, one could as well resolve it with a precision 
measurement of sin^ 2^13 with a reactor experiment. The JHF-SK setup is in Figure 10 
shown in three different configurations: alone with neutrino running only (JHF-SK), alone 
with the total running time split in order to obtain about equal numbers of neutrino and 
antineutrino events (JHF-SKec), and in combination with Reactor-ll with neutrino running 
only (JHF-SK-|-Reactor-ll). The combined neutrino-antineutrino running at JHF-SK^c clearly 
helps to improve the performance compared to running with neutrinos only (JHF-SK). There 
is, however, only a marginal improvement by adding the reactor experiment to JHF-SKcc- 
The reason is that the correlation between sin^ 2^^i3 and (5cp is in this case already resolved 
by the antineutrino channel, and a precision measurement of sin^ 2^^i3 does not improve the 
sensitivity to CP violation further. The third, most interesting case in Figure 10 shows the 
combination of JHF-SK with neutrino running only with Reactor-ll. Because the statistics 
for JHF-SK is much better than for the combined neutrino and antineutrino running here. 



22 



and the correlation between sin^ 2^13 and 5cp can be resolved by the reactor experiment, 
the overall performance is optimal. It even covers the LMA-II best-fit value, which leads 
to the interesting chance to observe leptonic CP violation by such a combination. We have 
not found any other combination {e.g., involving two superbeams as in Ref. [15]), which 
could achieve such a good coverage in the sin^ 2^i3-Am|i-plane by using the first-generation 
JHF-SK and NuMI superbeams. Especially, the CP performance of NuMI suffers from the 
fact that mainly the first term in Equation (1) is enhanced by matter effects, which means 
that the relative weight of the second and third CP-sensitive terms is lower than in the 
JHF-SK case. A reactor experiment could thus be a very important element to measure 6cp 
without including an antineutrino running at superbeams. This is, however, because of the 
limited statistics, for the first-generation superbeam experiments restricted to the LMA-II 
region. 

7 Summary and conclusions 

We have presented a detailed study of potential reactor experiments with near and far 
detectors. The first important aspect covered by this study has been a careful modeling 
and discussion of the systematical errors of such reactor experiments. The second major 
issue has been the comparison to the first-generation superbeams JHF-SK (JHF to Super- 
Kamiokande) and NuMI and the study of complementary effects between reactor experi- 
ments and superbeams. For a quantitative discussion, we have defined benchmark reactor 
experiments with far detectors up to the size of the KamLAND detector, which we have 
labeled as Reactor- 1 and Reactor-ll. They have integrated luminosities of 400tGWy and 
SOOOtGWy, respectively, where the units are given as detector mass [tons] x thermal 
reactor power [GW] x running time [years] . 

The structure of the relevant oscillation probabilities of reactor experiments demonstrates 
that they are only very little spoilt by correlations and not at all by degeneracies (c/.. Sec- 
tion 2). Very good sensitivity limits for sin^ 2^13 can therefore be obtained. For our setups, 
we find, for the current best-fit values, sensitivity limits of sin^26'i3 ~ 1.8 • 10^^ (Reactor-I) 
and sin^2^i3 ~ 7.6 • 10~^ (Reactor-ll), which is a significant improvement of the current 
CHOOZ bound. Compared to the superbeams, the Reactor-I sensitivity limit corresponds 
to the JHF-SK experiment including systematics, correlations, and degeneracies, which is 
about sin^ 2^^i3 ~ 1.7 ■ 10^^ for the LMA-I solution. The dependence of the sensitivity 
limits of Reactor-I, Reactor-ll, and JHF-SK on systematics, correlations, and degeneracies is 
illustrated in Figure 11. Compared to the superbeams, the sensitivities of the reactor exper- 
iments are almost unaffected by the true values of the solar parameters and less dependent 
on the value of Amg^. The dependence of the superbeam sensitivity limit on the LMA-I or 
LMA-II solution is shown in Figure 11, whereas the reactor experiments are hardly affected 
by Am|]^. The sensitivity limits of the superbeams are, for very large Am^i within the LMA- 
allowed region, not much better than the current CHOOZ bound. A reactor experiment is 
in contrast quite independent of the true values of the oscillation parameters. 

The careful modehng of the systematics leads to a number of important results. First of 
all, an appropriate choice of the energy cuts turns out to be very important. Especially, the 



23 



Sensitivity to sin^26i3 
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LMA-II 



10" 



Figure 11: The sensitivity to sin^2^i3 for Reactor-I, Reactor-II, JHF-SK at LMA-I (Aml^ = 
7 ■ 10"^ eV^;, and JHF-SK at LMA-II (Am\^ = 1.4 ■ IQ-^eV^; at the 90% confidence level. 
The sensitivitiy limits for the reactor experiments hardly depend on the true value of the 
solar parameters. The left edges of the bars correspond to the sensitivity limits from statis- 
tics only, the right edges of the bars to the real sensitivity limits after successively switch- 
ing on systematics (dark/blue), correlations (medium gray/green), and degeneracies (light 
gray/yellow). 



oscillation minimum should be within the energy window of the detectors, since spectral 
information is in reactor experiments very important for the reduction of systematical errors. 
The near and far detectors are assumed to be of identical technology in order to efficiently 
eliminate many systematical errors. It has turned out that the energy calibration error 
does not have a major impact on the results, since the spectral shape of the signal fixes the 
energy scale. The scaling of the sin^ 2^13 sensitivity limit with the integrated luminosity 
exhibits interesting features. In the low luminosity regime, we find that the value of the 
normalization error is rather important. Thus, the Reactor-I setup would lose approximately 
50% in sensitivity if the systematical errors are increased by a factor of 2. At large integrated 
luminosities, such as for the case of our Reactor-II setup, this dependence on the systematical 
uncertainty of the normalization is however much weaker. The reason for this is that 
with a high statistics sample it is possible to simultaneously measure sin^ 2^13 and the 
normalization with high precisions. Therefore, the sin^ 2^13 limit for the value of Unorm = 
0.8% is nearly the same as the one obtained with an unconstrained normalization, and, for 
luminosities above SOOOtGWy, the sensitivity limit again scales as At those large 

integrated luminosities, it becomes important if the calibration and normalization errors 
are sufficient to describe a real experiment, or if there are other sources of errors at the 
per mill level. We have therefore made a conservative estimate of the maximal admissible 
size of such additional error contributions and have found that an uncorrelated bin-to-bin 
error of 0.1% would not spoil the behavior up to a integrated luminosity of at least 
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10^ t GWy. Increasing this uncorrelatcd bin-to-bin error by a factor of five to 0.5%, however, 
would lead to a saturation of the sensitivity limit at luminosities around 5 000 1 GWy. For 
comparison, a detector of the KamLAND size at a 10 GW power station with a running time 
of 5 years would reach a luminosity ol C — 50 ODD t GWy. Those considerations indicate 
that a sensitivity to sin^ 2^13 of the order of 10~^ could be feasible with an ambitious new 
reactor experiment. It however remains an open question at which level the systematics of a 
real experiment will finally put a hard limit on the sensitivity reach of a reactor experiment. 
A sensitivity increase of one order of magnitude beyond the current CHOOZ bound down 
to sin^ 2^^i3 = 10~^ is however achievable quite independently of the size of the systematical 
errors. 

Even though reactor experiments are very good for sin^ 2^13 measurements, they cannot 
replace superbeams. Reactor experiments do not allow precision measurements of Amg^ 
independent of sin^ 2^^i3, and they are not at all sensitive to ^23- In addition, they do not 
have a significant sensitivity to the mass hierarchy and cannot access 5qv- Reactor exper- 
iments can, however, boost the performance of superbeams by the clean measurement of 
sin^2^i3, which helps to resolve correlations and degeneracies. One could use, for example, 
a large reactor experiment together with JHF-SK in the neutrino running mode only, in 
order to obtain a better sensitivity to leptonic CP violation than with JHF-SK in a com- 
bined neutrino-antineutrino running mode. The reactor experiment resolves in this case the 
correlation between sin^2^i3 and 5cp, and the higher event rates of the neutrino running 
mode improve the overall performance. Another example is the determination of the mass 
hierarchy. Neither JHF-SK nor NuMI alone could successfully determine the mass hierarchy. 
However, it has been shown in Ref. [15] that their combination with a longer NuMI baseline 
of 890km would be sensitive to the sign of ^m\^ at least for small values of /\m\^. It 
turns out that the Reactor-ll setup together with the two superbeams in this configuration 
could, furthermore, be sensitive to the mass hierarchy quite independent of the true value 
of Am|i, because Reactor-ll indirectly helps to resolve the sgn(Am|i)-degeneracy by mea- 
suring sin^2^i3 precisely. We do not observe this behavior with Reactor-ll in combination 
with only one of the superbeams. 

We conclude that the described reactor experiments are a very promising option on similar 
timescales to superbeams. Especially, the possibility to push the sin^ 26*13 sensitivity limit 
about one order of magnitude below the CHOOZ bound, down to sin^2^i3 ~ 10~^, is in- 
triguing. This may be important, since most neutrino mass models, such as texture models, 
are linear in flavor space, and it would be very surprising if the diagonalization predicted 
extremely tiny values of sin^2^i3. This means that the improvement of the CHOOZ bound 
achievable by reactor experiments would be very valuable to neutrino phenomenology and 
theory, and it would crucially influence the strategy of long-baseline discussion. For the 
superbeams, this implies that the task of finding sin^ 2^13 should be reconsidered, and that 
other observables, such as the mass hierarchy and CP violation, may get more weight in 
the planning. We have demonstrated that clever combinations of superbeams with a large 
reactor experiment may allow the determination of the mass hierarchy, and may even limit 
or measure (5cp with the first-generation superbeams. We therefore believe that the reactor 
discussion should be included in the optimization of superbeams. 
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A Systematical errors and the treatment of the near detector 



In this appendix, wc describe in detail the x^-analysis of the near-far detector complex and 
how we implement the systematical errors discussed in Section 4. 

We write for the theoretical prediction for the number of events in the ith energy bin of the 
near {A = N) and far {A = F) detector, respectively, 

T,^^(l + a + b^ + c,)N,^ + g^Mf, (8) 

and consider a x^-function including the full spectral information from both detectors: 



A=N,F 



Here, is the expected number of events in the ith energy bin of the corresponding detec- 
tor, which depends on the oscillation parameters, and Of is the observed number of events. 
We use 62 bins in the range between Ep^ — 1.8 MeV and Ep^ — 8.0 MeV, corresponding to a 
bin width of 0.1 MeV. In the absence of real data, we take as "observed number of events" 
the expected number of events for some fixed "true values" of the oscillation parameters. 
Per definition, the near detector is as close to the reactor as that no oscillations will occur, 
i.e., the Tf^ do not depend on the oscillation parameters, and we can set Of' = Nf^ . (We 
relax this assumption in Appendix C.) 

For each point in the space of oscillation parameters, the x^-function has to be minimized 
with respect to the parameters a, 6^, , g^, g^, and q modeling the systematical errors, 
such as described in Section 4.2 in items 1 to 4. The parameter a refers to the error on the 
overall normalization of the number of events common to both detectors, and is typically 
of the order of a few percent. Furthermore, the parameters and b^ parameterize the 
uncorrelated normalization uncertainties of the two detectors, where we assume that an 
error below 1% can be reached. The energy scale uncertainty in the two detectors is taken 
into account by the parameters g^ and g^. To this aim we replace in Nf- the visible energy 
-^vis by (1 -|- g^)Eyis- Then we have to first order in g^ 



Nf{g^) ^ Nf{g^ = 0) + / Mf with Mf - ""''^ 



dg^ 



(10) 



A typical value for this error on the energy calibration is a^ai ~ 0.5%. In order to model 
the uncorrelated uncertainty on the shape of the expected energy spectrum, we introduce 
in addition a parameter Cj for each energy bin. Note that all the parameters describing 
the systematical errors are at the percent level, which means that the linear approximation 
in Equations (8) and (10) is justified. Finally, the bin-to-bin uncorrelated experimental 
systematical error of item 5 in Section 4.2 is included with the help of the term containing 
(Texp in Equation (9). In this way we assume that the observed number of events in each 
bin and each detector Of has in addition to the statistical error ±1/ Of the (uncorrelated) 
systematical error ±ae^pOf. 
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In order to investigate the impact of the uncorrelated shape uncertainty (Jshapc and bin- 
to-bin experimental error (Xcxp, we have performed an analysis using the complete as 
given in Equation (9). However, as discussed in Section 4.2, our numerical results are very 
insensitive on the spectral shape uncertainty as long as cxshape ^ 2%. Therefore, it is save 
to set all of the c, to 0. Furthermore, we assume that (Jexp 0.1% can be reached, which 
means that setting this error to zero has little impact on the numerical results. Being left 
with the 5 parameters a, 6^, , , and , we can analytically minimize Equation (9) 
(with Ci = 0) with respect to the parameters and g^ . The minimum is located at the 
values 



b^'^aC,, g^'^aCg, (11) 



where 



' " " ^ {^ + M'al,){l+Afal)-M^alal,' ^ ^ 

with the abbreviations 

i i i ' 

Inserting Equation (11) into Equation (9), we obtain an effective x^-function for the far 

detector 

Finally, we can take into account that only the sum a + b^ appears in the theoretical 
predictions in Equation (8). Therefore, introducing the new parameter a — a + b^, we 
obtain 



with 



Of V^cal 

with 

Hence, we have an effective x^-function for the far detector which is given by Equation (17), 
where the information from the near detector is properly taken into account by the error 
on the normalization. The representative values (Tnorm = 0.8% and a^ai = 0.5% used in our 
calculations should be understood in the sense of Equations (17), (16) and (18). 

Let us eventually note that Equations (16) and (18) show the correct behavior in all limiting 
cases. For no near detector at all (A^/^, M^^ — > 0), we obtain C;,, — > and, as expected. 
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a — > cTq. For our calculations, wc assume the more interesting case of a very large number 
of events in the near detector A^cr? ^ 1. In this limit, ~ — 1, C„ ~ and 



^ ~ - + - • (19) 
< 



Assuming realistic values of Oa = 2% [7] for the flux uncertainty and Uf, = 0.6% 
for the detector-specific uncertainty, we obtain with Equations (18) and (19) an effective 
normalization error of (Jnorm — 0.8%, which is the value we have used for our numerical 
calculations. 



B Experimental details of reactor experiments and summary of 
key assumptions 

The detector technology and performance of our proposed reactor experiment is similar 
to the one of the CHOOZ [7] and KamLAND [4] detectors, and the Borexino counting 
test facihty [46]. Those detectors are based on a sphere filled with a liquid scintillator, 
which is separated by a plastic barrier from a buffer liquid between the sphere and the 
photo multiplier tubes. The typical energy resolution for such a detector is about (5 — 
10)%/V-^vis [4,47]. We are using an energy resolution of 5%/V-Evis [46], our results, however, 
do not change for 7.5%/\/S^, which is obtained for the KamLAND detector. In addition, 
we assume that the detector has a constant efficiency in the total analysis range from the 
energy threshold at Ey-^s — 1.0 MeV up to Ey^s — 7.2 MeV, which we divide into 62 bins 
corresponding to a bin width of 0.1 MeV. The normalization is chosen such that the event 
rate per unit thermal power of the reactor x fiducial detector mass x data taking time at 
a distance of 1 km is given by 227.5 events t~^ GW""*^ y~^ [47]. 

As shown in Ref. [47], the accidental background rate can be suppressed to one event 
per year for a 200 1 detector, which is completely negligible for our purposes. Therefore, 
we assume that it will be possible to construct a detector which is basically free from 
backgrounds. In addition, cosmic backgrounds can be efficiently rejected by using passive 
shielding of > 500 mwe (~ 200 m rock overburden) together with an active muon veto 
and pulse shape discrimination [47]. The number of events from geo-neutrinos is at most 
0.2 events t^^ y^^ [47], which gives a total contribution of at most 0.1% to the total rate 
and should therefore not present a limitation. For the short baselines used in this work, 
the contribution of other power stations is very small. Although some of the mentioned 
background sources may cause a total contribution of the order of 1%, their actual influence 
will be much smaller, since they are known to some extent and can also be measured during 
the periods the reactor is switched off. Thus, a 10% determination of a 1% background will 
contribute only to 0.1% to the total error, and neglecting the backgrounds is an excellent 
approximation for the far detector. 

In our standard setups we assume a near detector which is close enough to the reactor that 
no oscillations develop. (We relax this assumption in Appendix C.) In order to keep the 
characteristics of the near and far detectors as similar as possible it may be necessary to 
use exactly identical detectors. However, if it is possible to use a smaller near detector 
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it should have a size such that the event rate is at least ten times higher than for the 
far detector without oscillations. Of course, the near detector will have a much smaller 
rock overburden than the far detector, and therefore the cosmic background will be larger. 
However, because of the much larger event rate in the near detector the signal to cosmic 
background ratio should be even better than in the far detector, which means that we can 
also neglect the backgrounds in the near detector. Let us illustrate this by the following 
estimation: Consider, for example, our standard baseline of 1.7 km for the far detector, a 
baseline of 0.17km for the near detector, and a near detector size of a tenth of the far 
detector. For instance, with a rock overburden of 50 m ~ 125 mwe, the resulting muon flux 
would be higher by a factor of 10 [7]. Since the near detector has about a tenth of the 
volume of the far detector, its surface area is about 10^^^ ~ 4.6 smaller. Therefore, the 
total number of muon induced events is approximately twice as high as the one of the far 
detector. However, since the near detector has about ten times as many signal events as 
the far detector, the signal to cosmic background ratio should be a factor of 5 better than 
in the far detector, and neglecting the background is an excellent approximation. 

In the following we summarize the most important assumptions about the reactor neutrino 
experiment adopted in our calculations: 

• We consider one single reactor block. 

• We assume that neglecting backgrounds is a good approximation for the near as well 
as for the far detector (see Ref. [47] and the estimates above). 

• The full energy spectrum above the threshold is used. The impact of a low energy cut 
is investigated in Section 4.3. 

• Detection efficiencies are assumed to be constant in the full energy interval. 

• To convert the integrated luminosity L — fiducial detector mass [tons] x thermal 
reactor power [GW] x running time [years] into the number of events we follow the 
rule given in Ref. [47], assuming 

— a PXE-based scintillator, 

— the reactor running full time at nominal thermal power, 

— 100% detection efficiency. 

Other detector materials, a lower efficiency, reactor-off periods, or an operation at a 
lower thermal power, lead to a simple rescaling of our results. 

• We assume that the two detectors can be kept relatively calibrated during the full 
measurement period. 

• In our standard setups the near detector is situated as close (< 200 m) to the reactor 
as no oscillations develop. The impact of larger near detector baselines is investigated 
in Appendix C. 

• Finally, our assumptions about systematical errors and their relevance are summarized 
in Table 3. 
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Standard assumption 


Reactor- 1 


Reactor- II 


Effective normalization Cnorm = 0.8% 
Energy calibration (Jcai = 0.5% 
Exp. bin-to- bin uncorr. error cTexp ^ 0.1% 


important 
not important 
not important 


not important 

It 

not important 
important 



Table 3: Standard assumptions about systematical errors and their relevance for the two 
reactor benchmark setups Reactor-I and Reactor-ll used in this work. 

C The position of the near detector 

For practical reasons it might be hard to find a reactor station where a near detector can 
be situated very close (< 200 m) to the core with sufficient roclc overburden. Therefore, it is 
interesting to investigate the impact of larger near detector baselines on the sin^ 26*13 limit. 
In this case there will be already some effect of oscillations in the near detector and it is not 
possible to simplify the analysis as discussed in Appendix A. Therefore, we have performed 
an analysis based on the full x^-function given in Equation (9), taking into account the 
effect of oscillations in both detectors. Now the information provided by the near detector 
on the initial flux normalization and energy shape is already mixed with some oscillation 
signature. Hence, one expects that the correct treatment of the shape uncertainty cxshape 
due to the coefficients q in Equation (9) becomes more important. 



0.03 F 




Near detector position [km] 



Figure 12: The sensitivity to sin^ 29i-i at the 90% confidence level for Reactor-I and 
Reactor-ll as a function of the near detector position. The far detector is situated at 1.7 km 
and we assume identical detectors. Furthermore, the impact of an uncorrelated theoretical 
shape uncertainty of 2% is shown. 
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The results of this analysis are presented in Figure 12. We find that for the case of Reactor-I 
the hmit starts deteriorating around a near detector distance of 400 m, whereas for Reactor-I I 
the hmit even improves shghtly up to near detector basehnes of ~ 1 km. Due to the high 
statistics in the case of Reactor-I I, flux normahzation and shape are very weU determined by 
the near detector even in the presence of some effect of sin^ 2^13, and the additional infor- 
mation on oscillations improves the limit a bit. Furthermore, we find from Figure 12 that 
the shape uncertainty becomes important for near detector baselines > 1.1km, especially 
for Reactor-I I. A reduction of this theoretical error would be helpful in such a situation. 
We note that assuming the shape error to be completely uncorrelated corresponds to the 
worst case. A more realistic implementation of the shape uncertainty including correct cor- 
relations will lead to results somewhere in between the curves for no and 2% shape error in 
Figure 12. These calculations have been done for our standard value of Amg^ = 3- 10~^eV^. 
However, we find very similar behavior for other values of Am|]^. 

To summarize, for the case of Reactor- l-like experiments one should look for a site where 
the near detector can be placed at a distance of at most 400 m from the reactor. For large 
detectors, such as Reactor-I I, near detector baselines of up to 1 km will perform well. For 
near detector baselines longer than about 1 km the correct treatment of the theoretical 
shape uncertainty becomes important. 
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